Engineering solid-like structures via an arrested spinodal decomposition 
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The possibilities to tune the structure of solid like material resulting from arrested spinodal 
decomposition is investigated using a system composed of lysozyme, a globular protein, dispersed 
in a water solution as model system for colloids with short range attraction. It is shown that the 
resulting arrested spinodal decomposition is driven by the interplay between the early kinetics of the 
spinodal decomposition and the dynamical arrest. The initial concentration, the quench depth and 
speed from the fluid state to the arrested state enable to tailor the mesh size of the solid network 
of the arrested spinodal decomposition. 

PACS numbers: pacs 



I. INTRODUCTION 

Dynamical arrest reveals itself to be a powerful mech- 
anism to prepare solid-like material and is widely used 
in food science, care product industry and material sci- 
ence pQ. Recently, efforts have been made to understand 
the influence of interactions on dynamical arrest [21 [3]. 
In particular in short range attractive colloid-like system 
it is possible to form solid-like network through an in- 
terplay between dynamical arrest and phase separation 
[11 |5j [6j [7]. In such systems, the spinodal decomposi- 
tion process leads to a bi-continuous structure which gets 
'pinned' into a rigid self-supporting network because the 
particle-rich regions have undergone dynamical arrest be- 
fore reaching the equilibrium conditions set by the coex- 
istence curve. Experimentally this pinning of the phase 
separation due to dynamical arrest manifests itself as a 
crossover from a time dependent to a time independent 
behavior of the structure factor. Our ability to tune the 
properties of this arrested network depend very much 
on our understanding of its formation precisely because 
the resulting structure is out of equilibrium. We thus 
focus on the kinetic and path dependency (initial con- 
ditions such as the concentration in colloids and quench 
rate and depth from the fluid state to the arrested spin- 
odal decompositon) on the final structure of the arrested 
spinodal decomposition. 

In a first part we characterize the system, lysozyme, a 
globular protein, dispersed in water and develop the anal- 
ogy with colloids. In a coarse grain approach, scattering 
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data from the fluid phase, above the phase separation 
allows to extract interaction parameters of lysozyme in 
water such as the range and depth of the short range 
attraction. In a second part we show that for shallow 
quenches into the phase separation region, the system 
exhibits the classical features of spinodal decomposition 
whereas deeper quenches below the spinodal line lead to 
the formation of arrested spinodal decomposition where 
the resulting structure can be tailored mainly by con- 
trolling the rate and the depth of the quench. Finally 
we investigate the structure formed in the course of the 
arrested spinodal. We obtain quantitative and consistent 
information about the structure of the arrested spinodal 
with a combination of scattering and optical microscopy 
over a very large range of length scales. In particular, at 
low scattering vector moduli, q, (i.e. large length scales) 
we combine microscopy and ultra small-angle light scat- 
tering (USALS) to get a straightforward analysis of the 
correlation length and mesh size of the arrested network. 
Similarly, at high q (shorter length scales) we use small- 
angle neutron and x-ray scattering (SANS and SAXS, 
respectively) combined with a simple assumption of lo- 
cally liquid-like structure in order to obtain rather de- 
tailed structural information that is fully consistent with 
the position of the arrest line determined in [1]. However, 
the situation at intermediate q is much more difficult if we 
like to get information beyond the fact that an extended 
Porod region [8] with q~ 4 behaviour indicates the forma- 
tion of well separated regions of different concentrations 
with a well-defined interface typical of phase separation. 
Here we have introduced reverse Monte Carlo approach 
coupled to the coarse graining procedure to limit compu- 
tational effort in order to obtain additional information 
and describe the problems that go hand in hand with the 
very large range of length scales. 



2 



II. MATERIALS AND METHODS 

The study of the structures formed during the spin- 
odal decomposition process requires the use of a broad 
variety of experimental techniques to bridge the gap be- 
tween the building blocks, the lysozyme molecule, which 
sets the lower limit of the required accessible length scale 
to the nanometer, and the resulting final structure, a 
network with mesh sizes of a few microns which deter- 
mines the upper limit [5J [TU1 E] • The large scale struc- 
tures are probed by phase contrast microscopy or by ultra 
small angle light scattering. The lower scale organization 
is monitored using small angle neutron scattering, small 
angle x-ray scattering and static light scattering (SLS). 



were taken so that the experimental results were compa- 
rable, reproducible and artifact free. The illumination of 
the sample is kept constant to compare different images. 
To minimize the heat effects of the illumination, we illu- 
minate the sample only during the time required to take 
pictures. Moreover we attenuate the heating effect of the 
UV and IR part of the illumination spectrum by plac- 
ing a beaker of water between the light source and the 
sample. The microscope focus was set to image the mid- 
plane of the sample to avoid the influence from wetting 
effects at the interfaces between the lysozyme solution 
and the capillary. The objective aperture was opened to 
its maximum diameter to minimize the depth of the focal 
plan to a few microns. 



A. Sample preparation 

We use hen egg white lysozyme (Fluka, L7651) dis- 
persed in an aqueous buffer (20 mM Hepes) containing 
0.5 M sodium chloride. Lysozyme is a monodisperse glob- 
ular protein of a molecular mass of 14.4 kDa carrying a 
net charge of +8 electronic charges at pH= 7.8, [T2]. Ini- 
tially a suspension at <j) w 0-22 is prepared in pure buffer 
without added salt, and its pH is adjusted to 7.8 ±0.1 
with sodium hydroxide j4j [13]. Under these conditions 
the suspension is stable and is further used as stock. We 
then dilute it with a NaCl-containing buffer at pH= 7.8 
to a final NaCl concentration of 0.5 M. Particular care is 
taken to avoid partial phase separation upon mixing by 
pre-hcating both buffer and stock solution well above the 
coexistence curve for liquid-liquid phase separation (cf. 
Fig. 1). This procedure results in completely transpar- 
ent samples at room temperature with volume fractions 
ranging from <fi — 0.01 to 0.18, where <f> was obtained 
from the protein concentrations c measured by UV ab- 
sorption spectroscopy using (j) = c/p, where p — 1.351 
g/cm 3 is the protein partial specific density, |14j . To 
prepare samples at even higher volume fractions up to 
<\> Ki 0.34, we take advantage of the ability of the system 
to phase separate into a protein-rich and protein-poor 
phase. Typically, a sample at = 0.155 is quenched to 
a temperature 15°C<T< 18°C below its cloud point and 
centrifuged at 9000g for 10 minutes. Quasi-equilibrium 
is reached once the two phases are separated by a sharp 
meniscus and show only slight turbidity, [I]. The super- 
natant is then removed and the bottom dense phase is 
used for further experiments. 

B. Microscopy 

The optical microscopy is performed with a Leica DM- 
IRB in phase contrast mode. The lysozyme solution is in- 
jected in a flat capillary tube (width 1 mm, depth 50 /im) 
while it is in a state of a homogeneous fluid at 25°C. The 
sample in the capillary is embedded in a thermostated 
holder on the microscope stage and imaged. Precautions 



C. Static light scattering 

Light scattering experiments were performed with a 
commercial ALV/DLS/SLS-5000 monomode fibre com- 
pact goniometer system with an ALV-5000 fast corre- 
lator. We used NMR tubes as scattering cells (5 mm 
diameter). The data were corrected for background scat- 
tering (cell and solvent) and converted to absolute scat- 
tering intensity - the so called Rayleigh ratio - using 
toluene as reference standard. The Rayleigh ratio can be 
written as / = KcMPS where K is the contrast, c the 
lysozyme concentration, M the molar mass of lysozyme, 
P the form factor and S the structure factor. The refrac- 
tive indices necessary for calculating the contrast term, 
which is given by K = (2irnodn / dc) 2 /(TV^A 4 ) with no 
the index of refraction of the buffer, dn/dc the refractive 
index increment, and the wavelength, A, of the scattered 
light in vacuum, were determined for all solutions with an 
Abbe refractometer by interpolating the measurements 
of the index of refraction, done for three different wave- 
lengths of 435.1 nm, 546.1 nm, 579.1 nm, to the ALV- 
5000 laser wavelength, 514.5 nm. The resulting refractive 
index increment, dn/dc=0. 194 mL/g, was found to be 
temperature and salt independent within an error of 2%. 
Since lysozyme is much smaller than the wavelength of 
the ALV-5000 laser we are in the long wavelength limit. 
We are thus assuming that the scattered intensity is q- 
independent , P(q) =1, that S(q) ) is obtained directly 
from the experiment. 

The USALS setup is described elsewhere 15]. It covers 
a q range of 0.1 to 2 /jm" 1 . The sample is filled into 
square cells with a short optical path length of 10/xm 
to avoid multiple scattering and then transferred into a 
thermostated cell holder. 



D. SANS and SAXS 

SANS and SAXS are ideal tools for studying the struc- 
ture of materials in the mesoscopic size range between ~1 
and ^400 nm. SANS experiments were performed at the 
SANS I facility at the Swiss neutron source SINQ at the 
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Paul Scherrer Institut, Switzerland, and at the instru- 
ment Dll at the Institut Laue Langevin in Grenoble, 
France. We used 1 mm Hellma quartz cells and a ther- 
mostatically regulated sample holder. Combinations of 
different wavelengths, sample-to-detector distances and 
collimation lengths were used to cover a q range of 0.015 
to 2 nm _1 . In the SANS experiments, the contrast, 
K 7 depends on the scattering length density of water 
Pl =-5.60xl0 7 Ar 2 , on the scattering length density of 
lysozyme pw = 1.88xl0 6 A~ 2 , the partial specific vol- 
ume of lysozyme, v — 0.74 mL/g [T3] and the Avogadro 
number, Na- It is equal to K = v 2 (pl — Pw) 2 /Na — 
5.41xl0~ 8 m 2 .mol.g~ 2 . SAXS experiments were carried 
out in Department of Chemistry at University of Aarhus, 
Denmark with a pinhole camera (NanoSTAR, Bruker 
AXS) employing a rotating anode (Cu Ka) x-ray source, 
a thermostatically regulated sample chamber and a two- 
dimensional gas detector |16j . For the present experi- 
ments the sample to detector distance was 24 cm provid- 
ing q range of 0.2 to 8 nm _1 . 



E. Temperature quenches 

The phase diagram of the lysozyme solution we use 
(Fig. 1) exhibits the classical features of a colloidal sys- 
tem with short range attractions [17] • We observe a 
liquid-liquid coexistence curve that is metastable with 
respect to the solid-liquid phase boundary. The dynami- 
cal arrest line intercepts the coexistence curve leading in 
particular to arrested spinodal decomposition [U [18] . We 
investigated two types of temperature quenches differing 
by their cooling rate. For both types of quenches, the 
cuvette used in the experiment was filled with a fresh 
lysozyme dispersion around 25° C. As shown in Figure 
1, the sample is in an homogeneous fluid state provided 
we work fast enough so that crystals don't have time to 
develop (crystal start forming after one to two hours af- 
ter the sample preparation). The lysozyme solution is 
then quenched below the spinodal line to the desure final 
temperature, TV. In the 'fast' quench the sample is pre- 
quenched for a minute in an cthanol bath at TV. It is then 
quickly transferred to the thermostated cuvette holder of 
the experiment, also at TV. This procedure allows to ob- 
tain the fastest experimental quenches possible with this 
cell size. We estimate the time for the sample to go from 
25°C to Tf to be less than 30 s. 'Slow' quenches result 
from a simple cooling down of the sample from 25° C to 
Tf as the sample is inserted in the pre-thermostated sam- 
ple holder of the experiment. It then takes about 100 s 
to thermalize the sample to TV. 



III. LYSOZYME AS A MODEL COLLOID WITH 
A SHORT RANGE ATTRACTION 

We worked with lysozyme which is a globular protein 
interacting via a temperature dependent short-range at- 
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FIG. 1: State diagram of lysozyme, Hepes buffer with a pH 
of 7.8, [NaCl]=500 mM. 



traction and a combination of a hard core and a soft 
long range repulsion due to surface charges (around +8e 
at pH=7.8). We prepared the lysozyme sample at high 
ionic strength, [NaCl]=500 mM/mL, where the hard core 
repulsion and the short range attraction are dominating 
due to the screening of the charges on the proteins. Under 
these conditions the phase diagram, Fig. 1, is in agree- 
ment with simulations and experiments done on colloids 
interacting via a hard core repulsion and a short range at- 
traction: the gas-liquid curve is metastable with respect 
to the fluid-crystal boundary and the dynamical-arrest 
line intercepts the coexistence curve leading to an ar- 
rested spinodal decomposition [HUH]. To further exploit 
the colloid-lysozyme analogy we measure the structure 
factor of lysozyme in a wide range of concentrations and 
temperatures. A comparison of the measurements with 
the structure factors calculated for different model inter- 
actions such as a DLVO potential, a square well potential, 
etc... allows us to obtain additional insight into lysozyme 
interactions. 

Since 500 mM of NaCl screens almost totally the 
Coulomb repulsion, we decided to work with a square 
well interparticle pair potential, Usw{ r ) to model short 
range attraction for lysozyme: 

oo < r < 2a 
Usw{r) = {s 2a <r <2a.(l + A) 
r>2a.(l + A) 

where a is an effective hard sphere radius, and e and A are 
the depth and the range of the square well, respectively. 
The square well potentials have unphysical shapes but 
they have been widely used to study colloids with short 
range attractions since the physics depends only weakly 
on the shape of the potential, [19 . 

In order to determine the rest of the parameters of the 
potential, e and A, we subsequently fit a temperature and 
a concentration series in the fluid state above the coex- 
istence curve, Fig. 2. The series contains 5(0) obtained 
from light scattering and S(q) obtained from SANS and 
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SAXS. The effective structure factor S(q) is determined 
as the ratio between two measurements normalized with 
the particle concentrations: 



5(c,q) 



I(q,c).Cdilute 
I{q,Cdilute)-C 



(3.1) 



where I(q, c) is the intensity scattered by the solution of 
concentration c and /(g, Cdiiute) is the intensity of the di- 
lute sample, Cdiiute=7 mg/mL, used as the effective form 
factor. The inset in Fig. 2b shows I(q, CdUute) measured 
both using SANS (0.015 nm" 1 < q < 2 ran" 1 ) and SAXS 



(0.2 nm" 1 < q < 



The theoretical expression for the effective structure 
factor was obtained through numerical calculations using 
integral equation theory with a Percuss- Yevick closure re- 
lation for a model of polydisperse spheres and an interac- 
tion potential given by a square well potential, |20|, I21j. 
We use polydisperse spheres to account of the slightly 
elliptical shape of lysozyme. This approach is further 
applied in the structure factor calculation. 

We obtain a satisfying fit of the temperature depen- 
dence of the structure factor with a square well interpar- 
ticle interaction. For the concentration series at a given 
temperature the range and the depth of the well could be 
maintained constant, Fig. 2a: e/fcT=-3,52 for 20°C and 
-3.24 for 30° C, A = 2.22%. This simple model under- 
estimates both the position of the critical point as well 
as the values of 5(0) at high volume fractions. Those 
effects might come from the fact that we have taken into 
account the ellipsoidal shape of lysozyme by introducing 
polydispersity on spherical objects. 

Moreover, the study of the structure factor as a func- 
tion of temperature, T, gives a direct relation between the 
depth of the short range attraction, e and T. For a given 
concentration, c—70 mg/mL, S(q) could be fitted main- 
taining the range of the attraction constant (A=2.2%) as 
shown in Fig. 2c. This is consistent with the presence 
of an isosbestic point at q OSnm" 1 . Isosbestic points 
indeed appear in scattering experiments when adjusting 
the strength of attraction at constant range, [2"2l [2"3"1 12"4"] . 
The analysis of the data shown in Fig. 3b provides an es- 
timate of the temperature dependence of e as shown in 
Fig. 2c. e/kT scales linearly with temperature, which 
would suggest that the main temperature dependence of 
the interparticle interaction is simply through the tem- 
perature factor, while other specific temperature effects, 
such as changes in the hydration of lysozyme, if present, 
play a minor role. In Fig. 3d, we compare the second 
virial coefficient measured from SLS and the one calcu- 
lated based on (e,A). We notice a similar temperature 
evolution but the calculated second virial coefficients un- 
derestimate the attraction. 

This study is in agreement with recent work on 
lysozyme where an adhesive hard sphere model had been 
used. Rosenbaum and Zukoski as well as Piazza et al. 
[231 l2"rT] found a similar temperature dependence of the 
potential well depth. Similar conclusions on the weak 
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FIG. 2: (a) 5(0) for 20°C (□) and 30°C (O) as a function of 
the volume fraction. The lines show the calculation results, 
(b) S(q) for different temperatures (18, 20, 25, 30, 40° C) at 
c=70mg/mL. S(q — > 0) increases monotonously as tempera- 
ture decreases. The inset show the experimental form factor 
of lysozyme. The high q data is corrected for internal scatter- 
ing length density fluctuations. The lines show the calculation 
results, (c) Evolution of e obtained from the fit of S(q) on (b) 
as a function of temperature, (d) Comparison of the second 
virial coefficients obtained from SLS (□) and the one calcu- 
lated (O) from (e, A). 
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temperature dependence of the potential had also been 
made by Malfois et al. [27] . 



IV. 'CLASSICAL' SPINODAL 
DECOMPOSITION 

For shallow quenches below the spinodal curve but 
above the arrest tie line at 15° C as shown in Fig.l, 
phase separation proceeds classically via spinodal decom- 
position. After some time, we obtain two homogeneous 
phases, a gas like phase with volume fraction </>i, and 
a liquid phase with volume fraction 02, separated by a 
sharp interface. 

We proceed with a 'slow' quench and observe the 
spinodal decomposition under the microscope. At early 
times, the micrographs in Fig. 3a show the typical evo- 
lution of a bicontinuous network with a characteristic 
size which increases with time, indicating a coarsening of 
the structures formed during the spinodal decomposition. 
After about t ^100 s, the dense phase starts to sediment 
due to the density mismatch. Micrographs indeed show 
an increassingly blurred appearance when the focus is 
moved to the upper part of the capillary. From then on 
the focus is maintained in the lower part of the capillary 
tube. Finally, around i=1000 s, the dense phase wets the 
bottom of the capillary and spreads to totaly cover the 
bottom of the capillary tube. A quantitative evolution 
of the correlation length, £, is obtained by directly mea- 
suring in the micrographs or by Fourier transforming the 
micrographs. The evolution of £ is presented in Fig. 3b 
and reflects the three regimes previously described: (a) 
growth (t < 100 s), (b) sedimentation (100 s< t <10000 
s), and (c) wetting-spreading (i > 10000 s). 

In parallel USALS experiments were carried out. In 
contrast to microscopy, where a small portion of the sam- 
ple is imaged, USALS, due to the fact that a large scat- 
tering volume is probed, allows to obtain better statis- 
tics. The evolution of the spinodal decomposition pro- 
cess and the corresponding characteristic length scale is 
captured in the Fourier space as shown in Fig. 4a. A 
peak is observed in the scattering pattern at a scatter- 
ing vector q* that corresponds to a characterestic length, 
£ = 2n/q*. With time the peak increases in amplitude 
and moves towards lower scattering vectors, indicating a 
coarsening of the structures formed during the spinodal 
decomposition. At t ^100 s, the evolution is perturbed 
either by sedimentation of the dense phase or by the fi- 
nite size of the cell. The scattering patterns shown in 
Fig. 5a are typical for spinodal decomposition seen for 
example in binary mixtures, polymer blends or in classi- 
cal colloid polymer mixtures [2~5J |2~H] . This is illustrated 
in Fig. 5b, where the data are shown to follow universal 
scaling typical for spinodal decomposition processes. We 
see that the scattering follows a form given by I(q)q* t 
as proposed by Furukawa [30] and derived by Dhont for 
the spinodal decomposition of colloids in the initial and 
intermediate stage including hydrodynamic interactions 
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FIG. 3: (a) Micrographs showing the time evolution of the 
spinodal decomposition for a shallow quench at 17°C, 1°C 
below the coexistence curve with an initial concentration of 
00=0.15. (b) Related evolution of the characteristic length, 

e- 



[31 . We see that all the data obtained in the beginning 
of the intermediate regime collapse onto a single univer- 
sal scaling function. The relationship £ = 2n/q* between 
the peak position q* and the characteristic length scale 
£ now allows us to investigate the temporal evolution of 
the domain size. Commonly one distinguishes between 
at least three characteristic regimes during the spinodal 
decomposition process [15]: (?) early, (ii) intermediate 
and (Hi) late stage [2lj]. The early stage is described by 
linear Cahn-Hillard theory [3"2l I3"3"] and in essence pre- 
dicts an increase of the amplitude at constant q* until 
the peak spatial concentration fluctuations have reached 
the final coexisting concentrations and domain growth 
and coarsening sets in. In a simple picture coarsening is 
expected to occur according to a simple power law £ ~ t a , 
where the scaling exponent is a = 1/3 in the intermedi- 
ate (or diffusive) and a = 1 in the late (or flow) stage of 
the coarsening regime [Ml [35]. A detailed calculation for 
the intermediate stage including hydrodynamics reveals 
that the exponent a is between 0.2 and 1.1, depending 
upon the relative importance of hydrodynamic interac- 
tions [31] . Our measurements capture the intermediate 
diffusive regime. Before any onset of perturbation effects 
such as sedimentation and finite cell size effects the scat- 
tering is in agreement with theoretical predictions and 
other experiments, (see [28 and references therein). 
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FIG. 4: (a) Static light scattering patterns obtained using 
small-angle light scattering from a phase separating sample 
for a shallow quench at f6°C, f.5°C below the coexistence 
curve with an initial concentration of <^o=0.11. As time goes 
the peak moves towards the low q: 0s (□), lis (a), 17s (o), 
34s (V), 101s (o), 118s (<), 151s (>), 250s (o), (b) Dynamic 
scaling of the data shown in (a). Inset: related evolution of 
the characteristic length. 



V. ARRESTED SPINODAL DECOMPOSITION 

We have previously demonstrated that lysozyme so- 
lutions at high ionic strength offer an interesting path- 
way to gelation via an arrested spinodal decomposition as 
sketched in Fig.l. If a sample is prepared above the crit- 
ical temperature T c but below the crystal line it remains 
in a homogenous liquid state characterized by the value 
of the volume fraction <^o and the temperature until fi- 
nally crystallization sets in. If such a sample is quenched 
into the spinodal region to a final temperature below the 
arrest tie line at 15°C the sample not only undergoes 
spinodal decomposition into a bicontinuous network with 
protein rich and protein poor domains, it also exhibits a 
non-ergodic transition. The non ergodicity transition is 
the result of dynamical arrest in the protein rich domains. 
The resulting arrested spinodal decomposition depends 
on the history of the sample. The history is indeed de- 
fined by the initial equilibrium state, in particular the 
concentration of the fluid sample, the rate of the temper- 
ature quench, the final temperature and the age of the 
sample. The quench rate influences the time the spinodal 
decomposition can proceed before it is frozen by dynam- 
ical arrest. Since the volume fraction of the arrested pro- 
tein rich phase, 4>2. glass, is lower than the expected value 
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FIG. 5: a) A sketch of the time evolution of density fluctu- 
ations for a system undergoing spinodal decomposition. At 
a certain moment the volume fraction (j>2 of the dense phase 
intersect the dynamical arrest threshold (f>2,Giass so that the 
phase separation gets pinned into a space spanning gel net- 
work. The initial and equilibrium volume fractions of the 
concentrate and dilute phase are denoted 0o, <fe and <j>\. 

4>2, we assume that the spinodal decomposition gets ar- 
rested during the early stage (Fig. 5). In this theory, [55] . 
fluctuations of concentrations initially have an amplitude 
that increases with time without changing wavelengths. 
The wavelength initially favored, q* , can be deduced from 
the amplification factor, R(q) |36] which depends on the 
relative values of the derivative of the chemical poten- 
tial and the energy cost of a concentration gradient. £ is 
related to to q* by £ = 2n/q* . 

We want to address the question: can we show quali- 
tative agreement between the theory of the early stage of 
the spinodal decomposition and the relevant experimen- 
tal parameters? By this we mean: can we show agree- 
ment between the evolution the structure and thus the 
characteristic length of the arrested spinodal decompo- 
sition, £, as a function of the final temperature, X/, the 
initial concentration, 4>q, and the position of the dynam- 
ical arrest line, 4>2, g iass(T)7 

We first focus on the effect of the overall concentration 
4>o at Ty=10°C and on the effect of T/ at constant </>o, 
Fig. 6. The quantitative evolution of the the characteris- 
tic length for these two series shows two trends. In the 
concentration series at constant temperature the char- 
acteristic length is constant or decreases slightly with an 
increase in initial concentrations. This suggests that R(q) 
is malmost independent of </>o as the arrest condition is 
the same: </> 2 ,gz q;!S (10 o C)=0.32. 

In the temperature series the characteristic length de- 
creases with decreasing temperatures (Fig. 6c). It im- 
plies that the wavelength initially favored, q* , moves to 
higher q values as Tf decreases. We also observe that the 
amplitude of the peak (Fig. 6c), diminishes with decreas- 
ing temperatures. Fig. 7 shows similar results obtained 
in USALS. This implies that the contrast is lower and 
that the phase separation gets arrested earlier when Tf 
decreases. This is in agreement with the fact that the dif- 
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FIG. 7: a) Static ligth scattering intensity obtained using 
USALS at </>o=0. 11 for different deep quench temperatures 15 
(□), 13 (o), 10 (A) and 8°C (V). b) Dynamic scaling of the 
data shown in (a). Inset: related evolution of the character- 
istic length. 



FIG. 6: (a) Micrographs showing the temperature and the 
concentration series, (b) Power spectrum concentration se- 
ries related to the micrographs to the data shown in (a). The 
series are performed for deep quenches at 10° C with 4>o equal 
to 0.075 (A), 0.11 (V), 0.15 (o), 0.19 (<). (c) Power spec- 
trum temperature series related to the micrographs to the 
data shown in (a). The series are performed at 0o=O.15 for 
quench temperatures equal to 1 (a), 5 (o), 10 (A) and 15°C 
(o). Insets in figures (b) and (c) show the evolution of the 
characteristic length respectively as a function of concentra- 
tion and temperature, respectively. 
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ference between <p2.Giass and </>i, which set the contrast, 
decreases as Tf decreases. 

We then look at the influence of the rate of the 
quenches on the characteristic length of samples show- 
ing arrested spinodal decomposition (Fig. 8). We ob- 
serve larger characteristic lengths for slow quenches at 
constant (f>Q. This in agrement with the previous as- 
sumptions. Since q* is shifted toward larger q values as 
Tf decreases and the system spends more time at high 
Tf in slow quenches and possibly arrests before reaching 
the final temperature Tf, and we expect to obtain larger 
characteristic length for 'slow' quenches. 

To conclude, the structure of the system created in the 
course of an arrested spinodal decomposition is mainly 



FIG. 8: Comparison of the evolution of the characteristic 
lengthes obtained from microscopy and USALS experiments 
as a function of Tf for two type of deep quenches, fast and 
slow quench. 



determined by the interplay between the early stages of 
the demixing process and the position of the glass line. 
Moreover, the resulting characteristic mesh size of the 
network is strongly influenced by the quench rate. Theo- 
retical support would be needed to test the following as- 
sumption: R(q) predominantly depends on Tf but only 
slightly on 4>q. 
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VI. LOCAL STRUCTURE OF THE ARRESTED 
SPINODAL DECOMPOSITION 

A. Simple approach 

In a next step we take a closer look at the local struc- 
ture of the arrested spinodal network probed with SAXS 
and SANS. We aim at characterizing the interface be- 
tween the gas-like and the glass-like phases and the local 
structure of the glass-like phase. The local structure of 
the dense phase depends on the local volume fraction 
but also on the interactions between the proteins. Given 
the fact that we expect the local structure to exhibit 
only weak variations as the system dynamically arrests 
when crossing the arrest line, we use liquid state theory 
to calculate the static structure factor at length scales 
comparable to the monomer diameter. We then compare 
the scattering results with model calculations where we 
use a square-well potential as defined in section III to 
approximate the protein-protein interactions. 

Figure 9a shows the evolution of the intensity as a func- 
tion of the final temperature Tf after a fast quench at 
00=0.15. Figure 9b shows the evolution of the intensity 
as a function the initial volume fraction, O at 10° C. At 
low q we observe a Porod regime, / ~ q -4 , which con- 
firms the idea that the gel is a bicontinuous network with 
a sharp interface between the gas phase and the dense 
glass-like phase. As temperature decreases, the intensity 
of the Porod regime increases, reflecting an increase of 
the surface to volume ratio of the glassy backbone back- 
bone (see equation below) . At the same temperature but 
different concentrations the network seems to have the 
same characteristic length, which then appears to depend 
mainly on the common arrest condition, 0=0.32. One 
can deduce the surface to volume ratio from the Porod 
regime: 

2nAp 2 S 

Iporod(q) = g4 y (6.2) 

In this equaton S is the total interface area and V is the 
volume of the sample. Ap is the excess scattering length 
density and it is given by contrast between the dilute 
phase and the dense phase. The gas phase is character- 
ized by its volume fraction 0i and the relative volume of 
the gas network 1 — h (h is determined using centrifuga- 
tion experiments, j4j). The dense phase is characterized 
by its volume fraction 2 and the relative volume of the 
dense network, h. Ap can be approximated by 

Ap = [p L (f>2 + (1 - <i>2)Pw\ - [pl4>i + (1 ~ 4>i)pw] (6.3) 

The results of the calculation of the volume to surface 
ratio using equations 2 and 3 are tabulated in Tab.l. 
The results indicate that the ratio V/S increases as the 
initial concentration decreases. This seems reasonable as 
the characteristic length and the local concentrations are 
independent of the initial concentration. The ratio V/S 



TABLE I: Variation of ci, C2, h (obtained from the cen- 
trifugation experiments [4]), £ (obtained from the microscopy 
experiments shown in Fig. 7) and V/S, the volume to surface 
ratio (extracted from the Porod regime) as a function of Tf 
and Co. 
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10 
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10 


c [mg/mL] 


200 


200 


200 


150 


ci [mg/mL] 


34 


55 


65 


55 


c 2 [mg/mL] 


465 


444 


353 


444 


h 


0.33 0.37 0.47 0.24 


£ H 


3.5 


2.5 


1.9 


2.5 


V/S [nm] 


332 


183 


76 


326 



increases with temperature. This seems also reasonable 
as the characteristic length increases with temperature. 

At larger q the intensity reflects the local organization 
of the dense phase and the particle interactions. S(q) 
was again obtained solving numerically the Percus-Yevick 
(PY) theory with a square-well potential including poly- 
dispersity using the algorithm in |20l 121] . The square- 
well potential parameters were taken from the fits from 
the previous section based on the extrapolation of e(T) 
obtained in the fluid region to Tf. Assuming that the 
contribution of the gas phase (1) and the glass phase (2) 
are uncorrclatcd, the intensity was estimated as follows: 

I(q,<fa,T f ) = K<p MP{q)S cal (q,cf> ,T f ) (6.4) 

where the calculated structure factor, S ca i, is given by 
[37]: 

S C ai(q,(/)Q,Tf) = [ — Z^ 1 Si(q, (f>i,e)+ h p L S 2 (q, 2i e)] 

00 00 

(6-5) 

Table 1 displays the values of the parameters used in 
the calculations. 

This approach is motivated by the fact that in general 
repulsive and attractive glasses show typical fluid struc- 
ture as the system goes through the dynamical arrest 
transition. The results of these calculations are shown in 
Fig. 10. Although this empirical model provides surpris- 
ingly good agrement at intermediate and high q, devia- 
tions at large q show that the approach does not capture 
the structure on larger length scale. Therefore the anal- 
ysis in the next section is applied. 



B. Russian Doll Reverse Monte Carlo analysis of 
the intensity 

We have shown that we obtain quantitative and con- 
sistent information about the structure of the arrested 
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FIG. 9: a) Temperature series: 20 (□), 13 (o), 10(a), 5°C 
(V). b) Concentration series: 150mg/mL (o) and 200mg/mL 
(a). (+) show the intensity of a dilute sample, 7mg/mL, at 
20°C. It is the reference for an effective form factor. The 
SANS data corresponds to points between q—1.5 10~ 2 to 2 
nm -1 . The SAXS data corresponds to points between 0.2 to 
8 nnT 1 . 

spinodal with a combination of scattering and optical 
microscopy over a very large range of length scales. A 
remaining difficulty is the situation of structure at inter- 
mediate length scales or wave vectors, in particular if we 
like to get information beyond the fact that the extended 
Porod region indicates the formation of well-separated re- 
gions of different concentrations with a well-defined inter- 
face. Therefore, we have developed a new reverse Monte 
Carlo (RMC) method for the analysis of the structure 
of the concentrated phase based on the scattering data 
(Fig. 11). 

Use of standard direct Monte Carlo would be extremely 
time consuming, as typically some 10 6 particles make 
up the biggest structures. Indeed, the experimental q 
range is very large (down to 0.01 nm -1 ), and the size 
of the primary lysozyme quite small (i? — 1.6 nm). To 
limit the computational efforts, we need a course grain- 
ing procedure, which we term 'Russian Doll RMC The 
idea is to describe the structure on a small length scale 
with few particles, regroup them into a 'meta-particle', 
and use it to build a higher order structure, and so on. 




0.1 1 0.1 1 

q [nm ] q [nm ] 



FIG. 10: Structure factor of the arrested spinodal de- 
composition, a) Quench at 13°C with co=200mg/mL. b) 
Quench at 10°C with co=200mg/mL. c) Quench at 5°C with 
c =200mg/mL. d) Quench at 10°C with c =150mg/mL. Dot 
lines stand for Si, dash lines for S2 and lines for S ca i- 



Each level is described by an ensemble of particles, us- 
ing a conventional reverse Monte Carlo (RMC) algorithm 
[551 159"] . The idea is to build a first-guess-structure, and 
improve the agreement of its I(q) in the corresponding q- 
range with the experimentally measured one by randomly 
displacing individual subunits. These displacements arc 
confined inside the spherical next order metaparticles, 
which is a way to control internal volume fraction and 
interactions. Due to possible interpenetration of the less 
dense metaparticles, the excluded volume condition is 
maintained only on the smallest length scale. Structure 
factors can then be calculated using the Debye formula 
[40 . Intensities are obtained by straightforward multipli- 
cation of the lysozyme form factor and structure factors 
for (on average) spherically symmetric metaparticles, fol- 
lowed by addition of the dilute-phase intensity. The in- 
tensity is average over many simulations to overcome one 
important short coming of the approach: the calculated 
scattering curve is that of only one structure whereas in 
reality, an ensemble of structures contribute. 

We focus on the temperature quench to 13°C (co=200 
mg/mL), but a very similar approach could be applied 
to the other samples. In Figure 11, the experimental 
intensity is compared to the successive fits on the var- 
ious length scales. At the smallest scale, 50 lysozyme 
particles make up a first metaparticle of radius 11.5 nm. 
Some deviations can be observed, presumably due to non 
sphericity of the proteins. At intermediate q, 70 of these 
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FIG. 11: Intensity scattered by the arrested spinodal decom- 
position (□). The thin lines show the intensity scattered by 
the meta-particles at different level. The thick line shows the 
result of the simulation. The Inset shows an example of the 
simulation's results in real space: the aggregation of third 
level metaparticles (radius 131 nm). 

metaparticles form the next higher structure of radius 44 
nm, 50 of which are then regrouped on the third level 
(R = 131 nm), and finally 200 of these biggest particles 
make up the largest structure (1000 nm). This structure 
represents 3.5xl0 6 lysozyme particles which make up the 
dense phase (35% volume fraction) inside the dilute one. 
For illustration, this biggest structure is plotted as an 
inset in figure 11. It reproduces the microscopic phase 
separation, on the typical length scale of which - two to 
three metaparticle diameters which approaches one mi- 
cron (consistent with light scattering and microscopy). 
The scattered intensity is found to be well fitted by the 
RMC simulation. In particular, all major features (low- 
q Porod, break in slope, high-g lysozyme structure) are 
reproduced in an at least semi-quantitative manner. 

To summarize, we have shown that our new Russian 
Doll-RMC leads to semi-quantitative agreement and rea- 
sonable consistency with other experiments. For future 
work, the limits of this approach will need to be tested, 
and in particular the set of working parameters (size and 



number of metaparticles) will have to be determined. 



VII. CONCLUSION 

We have exploited the protein-colloid analogy to show 
that the short-range attraction of lysozyme at a high salt 
content could be successfully modeled by a square well 
potential. In such systems we have shown that it is possi- 
ble to tailor the mesh size of an arrested spinodal network 
with respect to the quench rate and the quench depth. 
The origin of this comes from the interplay between the 
early stage of spinodal decomposition and the position of 
the glass line. This opens up new routes to prepare gel- 
like systems with adjustable structural and mechanical 
properties. It should be particulary useful in materials 
and food science. Moreover, the local structure of the 
arrested spinodal decomposition could be qualitatively 
modeled using a square well potential and an empirical 
decomposition of the structure factor into a contribu- 
tion from the dilute gas-like and concentrated glass-like 
phase. We have developed a new method to analyze the 
intensity scattered by arrested spinodal decomposition: 
Russian Doll Reverse Monte Carlo. The first results are 
encouraging (semi-quantitative agreement and a reason- 
able consistency with other experiments) and we will in 
the future put more effort into improving this analysis 
method. 
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